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Abstract 

A methodology is given to test the QCD Nf—2 chiral transition, presently conjectured to be 
second order. Scaling forms for the correlation length, susceptibilities and equation of state are 
given which account for finite lattice spacing. Confirmation by lattice simulation would provide 
a large set of consistency checks for establishing that the transition is second order. Further 
corrections from finite volume effects and higher dimensional operator mixing are given. The 
implications of scaling corrections in finite temperature lattice QCD studies are examined with 
emphasis on tests for the believed second order chiral transition within the realistic setting of 
finite lattice spacing effects. 
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Section I. Introduction 

Landau-Ginzburg-Wilson(LGW)-Hamiltonians have played a central role in the order parame- 
ter, renormalization group approach to second order phase transitions. Their importance to QCD 
was realized by Pisarski and Wilczek |j|], who pointed out that the QCD chiral phase transition 
has a suitable order parameter from which an effective critical theory could be constructed. By 
following standard notions of universality, they accounted for all the symmetries of QCD in the 
chiral limit in order to build an order parameter based effective theory. This led them to the class 
of SU(n)xSU(n) cr-models. Then, according to the universality hypothesis, it follows that these 
effective theories would give the same scaling behavior as the true microscopic theory, here QCD, 
near the critical point, provided that their order parameter ansatz is correct. 

Recently Wilczek [Q] has examined the case of two fermion flavors in the chiral limit where the 
effective LGW-theory is now the SU(2)xSU(2)=0(4), Heisenberg model. Since the latter theory 
has been well studied in statistical mechanics, the necessary information regarding the critical 
behavior of this theory was immediately available. Subsequently Rajagopal and Wilczek |^] have 
given a transcription dictionary which relates the 0(4)-Heisenberg transition to the language of 
lattice gauge theory. This transcription readily allows testing the original assertion, that the chiral 
phase transition is second order, by checking the universal scaling behavior and comparing to 
theory. 

This paper provides a theoretical consideration of critical phenomena scaling in the context of 
lattice gauge QCD. Although the general theory of scaling in critical phenomena is well studied, 
several attributes specific to lattice gauge theory have not been adequately appreciated. It is our 



goal to provide some of the missing ingredients here. As one outcome of these considerations, 
we hope it will motivate more extensive numerical lattice studies near the transition or crossover 
region, so that reliable scaling studies can be made. For this, a reasonable sample of points would 
be needed within a region no greater than 10 percent of the believed transition temperature. For 
example for the magnetization, < -0^ >, where simulation data is fairly reliable, presently there 
are about three simulation points within the fall-off portion of the magnetization profile. With 
this data, at best one can approximate the transition point. If one wants to conduct a serious 
renormalization-group scaling study, one needs at least a factor of five increase in data points in 
order to make reliable curve fits. Irrespective of how motivated the reader is by this paper, one can 
also turn to condensed matter physics, where present wisdom considers scaling studies as among 
the optimum tests for second order behavior. A resolution of the QCD chiral transition could now 
also benefit by intensifying quantitative efforts towards this means. Some efforts along these lines 
recently have been made by Karsch Q. He has examined scaling behavior of what he refers to as 
the pseudo-critical coupling. 

Theoretically, the considerations of |^, ^ regarding scaling are not complete since they overlook 
lattice specific corrections to leading scaling. Intrinsically, the effective theory in [Q] represents 
only a truncated sector of the complete QCD Hilbert space. In particular it does not include the 
full effect of short distance quark and gluon degrees. Their presence renormalizes the parameters 
of the effective theory and induces irrelevant, higher dimensional, operator correction terms to it. 
This leads to intrinsic scaling corrections to the critical behavior of QCD. Furthermore, in a lattice 
simulation there are also corrections to leading scaling which are extrinsic to continuum QCD. These 



effects arise from the lattice spacing, which at leading order induces dimension four corrections. 
Thus if one accepts the arguments of |jl], §, ^], or similar such based on symmetry considerations, 
that the chiral transition is second order, its critical behavior from a lattice simulation will have 
corrections to leading scaling due to the presence of the lattice spacing. In the analysis of g these 
effects were not treated. 

In total a lattice simulation involves four scales which in ascending order are the lattice spacing, 
the intrinsic cutoff Aqcd, the correlation length and the lattice volume. Ideally, the intrinsic length 
should be much bigger than the lattice spacing. In this case corrections associated with critical 
phenomena would be independent of the lattice spacing and only the continuum properties of QCD 
would be measured. In practice, typical hadronic lengths are only a few lattice spacings in size and 
are strongly dependent on the lattice size. Thus, it is important to understand the effects of finite 
lattice spacing since they reflect upon lattice artifacts that one is forced to live with in any real 
simulation. 

Following in the spirit of universality, we will discuss in this paper those aspects of nonideal 
lattice complications that induce universal effects. Since such effects will invariably appear in real 
simulations, it is useful to understand what we can learn from the well developed renormalization 
group theory of scaling and scaling corrections. What we will see is that this approach provides a 
way to study the critical properties near the chiral transition with rigorous account for the lattice 
spacing. As such it opens the possibility for lattice studies to circumvent their inevitable problems 
with lattice spacing effects and obtain further answers regarding the chiral transition. In addition 
to offering a strong set of consistency checks regarding the transition, since lattice effects can be 



explicitly treated, the methodology is also appropriate for further exploration of this transition 
from data obtained by simulations. 

The paper is organized as follows. In section two, as a brief refresher, we will review the 
renormalization group theory of second order phase transitions. In section three, scaling forms 
are presented which include lattice spacing correction effects. Other corrections arising from finite 
volume effects and higher dimension operators are also given for completeness. Finally, we give our 
conclusions and discuss various applications where scaling corrections will be of use in analysis of 
simulation data. 

Section 2 

In this section we will review the application of the renormalization group to studies of infra- 
red fixed points. Imagine one is given a Hamiltonian which, due to universality, can be any of a 
number of choices so long as it satisfies all the symmetry properties and has the correct spatial 
dimensionality as that of the original system of interest. In this respect, for our problem of the 
chiral phase transition a suitable such Hamiltonian as stated above is the 0(4) Heisenberg model. 
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Recall that since this describes a statistical mechanical system, implicitly there is a lattice cutoff 
a = 1/A. We wrote the theory in continuum notation to keep the closest analogy to field theory. 
However when doing momentum space calculations, recall that all integrals have an ultraviolet 
cutoff. Observe also that in statistical mechanics the bare Hamiltonian is the actual physical 
theory in contrast to quantum field theory where the renormalized theory is the physical one. The 



goal now is to study the universal critical behavior of this theory using the tools of field theory. 

The philosophy of the field theoretic approach to critical phenomena is the same as the Wilson- 
Kadanoff blocking approach in that one wants to coarse grain the original theory and generate a set 
of effective Hamiltonians, each describing larger cells of the original theory by way of thinning the 
dynamical degrees of freedom. To implement this idea field theoretically, one first rescales spatial 
and spin coordinates to appropriate magnitudes. In order to fix ideas, we will view the rescaling 
as a redefinition of the coordinate system with the physical system defined by the initial prime 
coordinates. Thus a momenta p for example in the system scaled by A as in (2-2) below becomes 
p/A in the physical system. From this point of view, as one examines a correlation function, for 
example, at a fixed momentum p but in a succession of scaled coordinate systems, meaning A -^ oo, 
one is studying the p ^ behavior of the correlation function for the original system. To define 
the scaled coordinate system x in terms of the original (physical) prime coordinate system x', for 
the spatial coordinates, since typical distances near the critical point are large, one defines 

x' = Ax (2 - 2) 

where A ^ 1. By this transformation, the relevant magnitude of x will now be of order one. Based 
on engineering dimensions, near the critical point the field variables are expected to be small in 
magnitude. As such one defines, 

*1M = ^^ (2-3) 

One now substitutes these quantities into (2-1) to obtain 

H = Jd'xkj:[v<Pa-V<Pa + mo<l>l]+^(J:A I (2-4) 
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where the bare mass rriQ is related to the temperature by rriQ ^ T — Tq with Tq being the mean 
field transition temperature. The first stage of the coarse graining is realized here in the sense that 
distances of order one in the unprimed coordinates are large in the original primed system. At 
an intuitive level one can conceptualize the effect of this coarse graining as saying that the energy 
density sampled in a cell dx in the scaled system is A times greater than in the original system. 
In the context of the Hamiltonian, the result of this first stage of coarse graining can alternatively 
be viewed as that of redefining the bare parameters of the unprimed system with respect to the 
original primed system. In particular the four body interaction coupling is, 

\/ A4-d 

u, = ^^ = goK'-'' (2-5) 

so becomes infinite as A ^ oo for d < 4 (recall that the relevant physical dimension of interest for 
finite temperature QCD is d=3). Note here that the transformation acts only on the spatial and 
field coordinates and not on the parameters which in particular include the lattice spacing. Thus if 
one examines a correlation function at momentum p in the unprimed coordinate system scaled by 
A with lattice cutoff a, it is equivalent to the correlation function in the original primed coordinate 
system at momentum p/ K but with the same lattice spacing a. In practical terms this means one 
can compute the renormalization properties of correlation functions in the limit A ^ oo, as one 
is used to doing in field theory, and at the same time be studying the infrared properties of the 
correlation functions of interest. 

To study critical phenomena with this formalism, one must demand that the theory be inde- 
pendent of A in the limit A ^ oo. Recall that in statistical mechanics there is no physical reason 
for demanding the above condition nor does it a priori lead to a scale invariant solution. It should 



be regarded as an ansatz for obtaining potentially scale invariant solutions. In this respect, this 
field theoretic approach cannot predict the critical temperature nor even determine a priori that 
a second order phase transition exists. The most it can do is assume the existence of one and 
check a posteriori. The latter is achieved if an infrared fixed point is found which will be discussed 
further below. Given that a fixed point is found, then the power of the field theoretic approach is 
its greater accessibility to analytic treatment compared to the Wilson-Kadanoff blocking approach. 
The quantities that can be computed are the scaling behaviors of all correlation functions, which 
in particular are seen explicitly to satisfy known scaling laws. From this critical exponents can be 
obtained, as well as the equation of state. 

How one defines the renormalization of the theory is not specified nor unique. In this respect, if 
a particular renormalization procedure does not yield an infrared fixed point in perturbation theory, 
it does not mean the initial theory has no second order phase transition. All one can conclude is 
that the sought for fixed point may require a different renormalization procedure in order for it to 
be reached by perturbation theory or that it may only be reached by nonperturbative means. On 
the other hand, if a fixed point is found, one has reasonable confidence that a second order phase 
transition exists in the system. 

Our discussion of renormalization will follow the conventional approach. Here one aims to take 
the A ^ cxD limit while holding fixed a certain minimal set of computed quantities. The statement 
of renormalizability in this standard approach is that one can adjust a finite number of parameters 
to make all observables finite. To renormalize the theory, we define the renormalized correlation 
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functions as, 

T(''Hp,go,A) = Z-^/Hgo,-)r^^\p,g,f^,A) (2-6) 

where Z are the renormahzation scahng parameters. The next step is to give the renormahzation 
conditions which we take as, 

^rg^ {p, g, fi\p2=^2 = 1 (2-7) 

^R\pi,g,fj')\sp = fj-^g 

The A ^ cx) limit can now be taken, and the renormahzed correlation functions will by construction 
be independent of A. Note that in implementing the renormahzation, a dimensionful parameter n 
is introduced. This is needed in order to specify the scale at which the renormahzation is done. 
The upper limit on fi is fixed by the lattice spacing, so that the study of lattice spacing corrections 
can be done with respect to a finite ^ as well as A. 

At this stage one must recall that our goal is to obtain the universal scaling relations. As such 
it does not matter what the specific initial bare parameters are. One is also free to analyze the 
theory in terms of bare or renormahzed quantities. To obtain the scaling relations the idea is, 
having fixed the bare parameter at cutofi' A, in order to maintain the renormahzation conditions 
(2-7) these parameters will develop a dependence on A. To determine their behaviors one uses 
the fact that the renormahzed theory is independent of A when A -^ cxd. Equivalently, one can 
demand that the bare theory be independent of fi. The role of /j, is similar to A in that ^u ^ oo 
corresponds to the limit of zero lattice spacing. Although the direct physical interpretation in terms 
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of renormalized quantities is obscured, it affords greater calculational ease and accuracy. Bearing 
in mind that even the bare theory is only an effective coarse grained version of the original lattice 
theory, not much is lost by this transcription. However if one prefers one can always keep in mind 
when making physical connections that the renormalized parameters are related by equations (2-7), 
or in different schemes by similar such, to the bare parameters of the actual system. 

Turning to the renormalization group equations, we will only discuss them to the extent of 
making precise our physical discussion above. Keeping in the spirit of this conceptual presentation, 
our discussion will be in reference to the bare theory, where physical contact is more transparent. 
For a more detailed account, the reader is referred to the very careful explanations given in |5|, |^. 
The demand that the renormalized theory be independent of the cutoff A gives, 

r('^)(p„5o,A) = (2-8) 



\ d ^. A. d N . A, 

oA /i ago 2 fi 



where 



Ad 
P{go,e,-) =A—\g^f,go (2 -9a) 

7]igo, e, -) = -A— |g,>Z i^go, -j (2 - 96) 

The solution to which is, 

r»(..,A,^.(..A-f^)e.p^f*l;f^ (--) 

The solution has the property 

rW(p,5o,^A) = r('^)(p,9o(^),A)expf r%(go(7)) (2-11) 

2 Ji 7 

which implies the A — > oo behavior is the same as studying the flow of ^(k;) as k ^ oo. To study 

the flow of g requires the study of the beta-function(2-9a). If (3 has a zero, the asymptotic limit of 
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g as K, ^ oo will be this zero provided f3'{g*) > 0. For a general coupling g ^ g* the solution (2-10) 
has nontrivial dependence on the scale parameter. On the other hand, at a fixed point, using (2-11) 
one obtains, 

r(^)(^,5o, A).^oo - «-'^+^('^-2+'')rW(p,5o*, A) (2 - 12) 

thus attaining scale invariance. 

Only at the fixed point is the theory scale invariant. To reach the fixed point requires either 
choosing the bare parameter to be exactly gQ or by taking the lattice spacing to zero and moving 
along the renornialization fiow trajectory. The former is highly unlikely whereas the latter is 
physically only approximately possible since statistical mechanical systems have a lattice cutoff. As 
such, corrections to scaling will occur and some parts of it are universal. In particular universality 
only applies to quantities computed at the fixed point. On the other hand, the imposition of 
a lattice, due to the finite lattice spacing, implies this fixed point is never reached under flow. 
Deviations from exact scaling are therefore reflected by deviations of the coupling from its fixed 
point value. For small deviations an expansion of the correlation functions about the fixed point 
with respect to the coupling constant can be done. From this Taylor expansion, only the derivatives 
of correlation functions evaluated at the fixed point are universal. This is the basic idea used to 
compute the desired corrections. In the next section we quantify this discussion. 

Section 3 

Our discussion of critical behavior will be in the context of lattice gauge theory. For this let 
us first relate the 0(4) theory to the effective QCD chiral Lagrangian as hypothesized by [|l|. The 
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order parameters of the effective chiral Lagrangian were stated to be the quark condensates, which 
can conveniently be written in matrix form as, 

Mi^i^qR-) (3-1) 

The effective Lagrangian must be symmetric to SU{2)l x SU{2)ji x U{1)l+r transformations. This 
can be expressed by parametrizing M with (a, it) and using the Pauh matrices to get 

M = a + iTr-f (3-2) 

Using this form for M in the effective Lagrangian, one can identify four real independent fields 
with respect to which the theory is 0(4) symmetric. 

Among the thermodynamic properties of the theory, the susceptibilities for T < Tc are identified 
with the masses of the a-particle and pion. There are two susceptibilities since symmetry breaks 
only in one direction. Defining a as the symmetry breaking direction, m^ is the longitudinal 

susceptibility, 

o dH , , 

m. = ^ (3-3a) 

and rn^ is the transverse susceptibility, 

^1 = ^ ■ (3 - 3fe) 

Let us now turn to the scaling behavior of the thermodynamic quantities and their correction 
from scaling. In the subsections to follow, we will discuss the three types of scaling corrections 
that will occur. These will be finite lattice corrections in subsection 3A, finite volume effects in 
3B and corrections from nonleading operator insertions in 3C. The former two give well prescribed 
effects that are not special to the particular choice of lattice action. An effective use of the latter 
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requires more specialized knowledge of the lattice action. The predictions are still in the context 
of universal scaling effects. However the usefulness of an analysis of operator corrections relies on 
first specifying an appropriate, manageable selection of terms that dominate. 

3A. Finite Lattice Spacing 

To obtain corrections due to finite lattice spacing, as discussed in Section 2, we must evaluate 
the solutions of the renormalization group equations at g away from g* . The basic idea is to Taylor 
expand the correlation functions about the fixed point and retain only the linear order deviation. 
This introduces one new universal exponent, lo = P'{g*), first discovered by Wegner Q. Details 
specific to our results below are given in appendix A. 

With corrections included, one obtains the following scaling behavior for the thermodynamic 
quantities ||8|, ^. Defining t = \T — Tc\, for T > Tc the correlation length behaves as, 

^ = A^{a)t''' {1 + C^Lia)^'') (3-4) 

the susceptibility as, 

X = A^{a)t-^l + C^L{a)tn (3-5) 

and the two point correlation function as, 

r(2)(p,5,t,/x) = ^2(a)/-'' [l + C2L(a)r" + Pip)] (3 - 6) 

For the equation of state, H{M) ,when T < Tc, there are two regions to consider, M/t^'^ <^ 1 
and M/t^'^ ^ 1. In the first case we have 

H = AH{a)M'f I ^ I (1 + CuLiay) , (3-7) 
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where the Widom function f{x) is defined in ^^, so that the longitudinal and transverse suscep- 
tibilities are respectively 

ml = A^ia)ml^{tr,Mr)il + Cj^L(a)r'^) (3 - 8a) 

ml = Arr,{a)ml(tr,Mr){l + CHLia)^'') (3 - 86) 



where m^j^o and mo-o are the leading scaling forms also given in [|T0|. In the latter case when 
M/t^'^ S> 1, the equation of state is, 

H = A^{a)M^ (l + CfjL{a)M'f^ (3-9) 

The r subscript on t and M in eqs. (3-7) and (3-8) is to signify the "rescaled" temperature and 
magnetization, which are defined in appendix A. These shift the absolute magnitude of t and M at 
different lattice spacings but not the scaling behavior with respect to them at a fixed lattice spacing. 
In testing lattice corrections to scaling, since the absolute magnitude of any one of the quantities 
is generally fixed at each lattice spacing by the leading scaling term, the precise knowledge of the 
reduced quantities is unneeded. Note that the rescaling factors affect all the temperature and 
magnitization variables above. However, where possible we have absorbed these factors into other 
coefficients, so as to leave the expressions in their most convenient forms. In the above the value 
of the critical exponents are u = 0.73 ± 0.02,r? = 0.03 ± 0.01, (3 = 0.38 ± 0.01, ,6 = 4.82 ± 0.05, 



7 = 1.44 lb 0.04, and w ~ 0.46. All but the latter are from a seven loop calculation in |13| and the 



latter is computed from a large-n expansion in p, 11, 12]. 

In the relations (3-4) to (3-9) the lattice dependence enters the correction term as the factor 
{g — g*) = L{a)t^'^ where we have explicitly represented the dependence on the lattice constant, 
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a, through the nonuniversal function L(a). Recall that the coupling constant g is an artifact of 
the specific renornialization scheme one adopts with no observable consequence. However since it 
enters in the same form, through the factor (g — g*) in the correction term, this implies the a- 
dependence of the correction term is reduced to the study of a single function. This is of relevance 
since one can empirically determine L{a) from a single measurement and apply the so determined 
result to all other quantities. For example it is possible to obtain a fit for the function L{a) from 
a given lattice measurement at two different temperatures of one of the thermodynamic quantities 
above. Then by the predictions of the theory this determines the function L[a) which controls 
the lattice dependence in the leading correction term in all the other thermodynamic functions. 
An alternative approach is to vary the temperature t, at fixed lattice spacing and fit the measured 
quantities to their respective forms in eqs.(3-4) to (3-9). One should then compute the ratio between 
the coefficient of f^^ to that of the leading term, for a given thermodynamic quantity, i, divided by 
the same for any other, j, which by the notation above is Ci/Cj. This should be the same at any 



lattice spacing a. In fact it has been shown that these ratios are universal |14]. Either approach 
would determine whether the predicted scaling corrections are there and thus further test for the 
existence of a second order phase transition. 

The definition of temperature T, in a lattice simulation requires some careful consideration. 
Recall that to vary temperature in a lattice with a fixed number of time slices, one varies the 
temporal lattice spacing. This in turn is accomplished by varing the coupling constant j3 = 6/5^ 
which is related to the lattice spacing through the /3-function. In such a procedure, therefore, not 
only does the temporal but also spatial lattice spacing change. If one is sufficiently within the 
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continuum limit, the intrinsic short distance scale would be independent of the lattice spacing and 
/?, or the temperature associated with it, could both serve as acceptable temperature parameters. 
This situation changes greatly in the more realistic case when one has not reached the continuum 
limit. Now as /3 is varied, both the temperature and correction coefficient, L(a), change. One 
remedy to this problem is to use different parameters in the temporal and spatial directions. If one 
is interested in studying detailed properties of continuum QCD, this is the most reasonable avenue 
to follow. However, there is a way to considerably simplify the task if all one is interested in is 
to establish the second order nature of the transition. If one accepts the assertion that continuum 
QCD with two light quarks is in the same universality class as the lattice-gauge action, then one is 
free to study the critical behavior of the latter in its own right without asking questions about how 
its parameters relate to QCD. In the context of critical phenomena, this follows since as far as the 
effective theory of the two systems is concerned, they differ at leading order in operator corrections 
by only differences in their respective short distance cutoffs. Thus, although the nonuniversal 
coefficients will differ, the leading as well as nonleading universal scaling behavior with respect to 
temperature is the same. With this simplification, one is free to use (3 as the temperature parameter 
of the lattice theory. Due to universality, this provides a pragmatic method for verifying scaling in 
the continuum theory also. 

Observe that the scaling corrections vanish at the critical point, t = 0. This is consistent 
with general notions about the second order phase transition where at the critical point there is 
only one length scale, that of the diverging correlation length, and the system forgets about all its 
microscopic length scales. Observe also that away from the critical point, t 7^ 0, any real many-body 
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system will exhibit scaling corrections since there is always a minimum microscopic length. 

Within the theory of [||, the two-point correlation function in (3-6) is supposed to be for the 
TT and a fields. Recall that in the context of the lattice theory, these fields were constructed from 
the fermion pseudoscalar and scalar bilinears respectively. In the theory of Q, the relation (3- 
6) is supposed to reflect the coarse grained, long distance behavior of these fundamental lattice 
correlators. Analogously, one can also write the two-point function for the p and ai fields which 
from [^ are, 

P\ = eimn(pmV4>n (3 - 10a) 

iA,)i = <A4V>z - </'zV>4 , (3 - 106) 

The two-point correlation function for these fields is defined as, 



{O,-''{x)Of\0)) ^ {6'^ - i|-)r(0'0'2)5„^,^<5 (3 - 11) 



with Oj"^ = e"^''' (/)'*' Vi0, 1 < a,/3, 7, 5 < 4, and has the scaling form with corrections, 

r(0'0'2)(^^i^^) = cok''-^g{kt-n{l + coo2L{a)rn . (3 - 12) 

The derivation of the correction factor also follows from what is given in appendix A. For these 
fields, accounting for corrections to scaling will be more important than for the vr and a fields since 
as pointed out by ^, mass terms dominate the nonanalytic part except very close to the transition 
point. In practical terms this means scaling will be difficult to detect so that maximal accuracy in 
the scaling form is desirable. 
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3B. Finite Volume 

Since all practical lattice studies are done on systems of finite size, studying the scaling behavior 
with respect to lattice size can be a useful tool. The renormalization group theory of finite size 
scaling [^, 15, 16] states that for any multiplicatively renormalizable thermodynamical quantity, 



P(t), such as the correlation length or susceptibility, 

^ = 5(i,i) (3-13) 

where 

g{L,t)=F{L/m) , (3-14) 

Phit) is the value of the thermodynamic quality in the system of finite size, and ^{t) is the correlation 
length of the same system in the infinite volume limit. Here F{x) behaves as, 

F{x) '^ 1 as X —> CO (3 — 15a) 

F{x) ~ x^/" as X ^ . (3 - 156) 

where p is the critical exponent associated with P{t). 

In a study which aims to test scaling behavior, the above can be used to compare two systems 
at sizes Li,L2, without having to know the L — > oo behavior. This can be implemented as follows. 
First measure, near the critical point, a quantity, PL-{t), as for example the susceptibility m{t), 
in systems one and two as a function of the temperature. From this and using (3-15b), find the 
associated temperatures fi,i2 such that 



L2J PL2{t2 
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If the transition is second order, then for any other multipHcatively renormaUzable thermodynamic 
quantity, the same relation (3-16) should be found at the same temperatures, except with p replaced 
by the critical exponent associated with the particular thermodynamic quantity under examination. 
Note that by this procedure, one never needs to determine the precise critical temperature in either 
of the two systems. In eqs. (3-16), for notational consistency, we used reduced temperatures on 
the right hand side. However, this relation can be applied with knowledge only of the absolute 
temperatures. 

3C. Operator Corrections 

The quartic 0(4) theory (2-4) is only the limiting coarse grained Hamiltonian obtained from 
the original lattice theory. This limiting Hamiltonian retains the most relevant operators and is 
exact for correlations in the zero momentum limit. Away from this limit, additive correction in 
an expansion in (p/A) arise for correlation functions at momentum p. The explicit form of the 
corrections has both universal and nonuniversal components. The universal part arises from the 
scaling behavior of the operator corrections. The nonuniversal part is the strength of each operator 
correction arising in the effective theory. The effective theory we are talking about here is the one 
that interpolates between the exact lattice action and the quartic LGW-Hamiltonian in (2-4). We 
write the effective Hamiltonian as, 

^eff — Hlgw + He (3-17) 

where Hlgw is (2-4) and contains the most relevant operators. He has all the less relevant operators 
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which we write as, 

oo Nj 

Hc = Y.ll<Ol (3-18) 

where j is the canonical dimension of Ol and the sum over i contains all operators in each set j. 
For the two-point function , operator corrections lead to the result, 

r(2) (p,g*, A) = r(2)(p,5*, A) + 5r(2)(p,5*, A) (3 - 19) 

where 

<5rW(p,5o*,A) =P^ (x ) '^"" (a)"'^ (^-'°) 

The theory of operator corrections has been applied for a variety of different purposes in studies 
of critical phenomena |^, g, 17, |l^]- For its relevance to lattice studies at finite temperature, both 



lattice specific and continuum QCD effects near the critical point may be further understood by this 
procedure. Such is an ambitious program which for one thing would need outside model specific 
considerations that limit the number of such correction terms to a small set, amenable to numerical 
fitting. By this decomposition one is able to systematically include operator corrections based on 
their order of importance. 

The dominant correction terms in H^ are from the operators of canonical dimension six. There 
are also dimension four operator corrections but their effect has been accounted for in the finite 
lattice spacing effects of subsection 3A. For the dimension six operators we have. 
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The anomalous dimensions they lead to are p^, 



r]i = m = m = 0, m = „ e, % = e, % = e, ??7 = 3e (C6) 



The theory of operator corrections suggests a flexible method for the analysis of critical behavior 
from lattice gauge simulations. In particular it will be of use if more extensive analysis of the fermion 
determinant is necessary. Our discussion in this subsection was meant only to be illustrative and 
to state its potential relevance for lattice gauge studies. 

Conclusion 

Model analogies and lattice studies have given plausible evidence that the chiral transition is 
second order, but as yet no proof. To make progress, increased quantitative efforts must be taken. 
As we have seen in this paper, the renormalization- group theory of critical phenomena provides 
an ideal possibility. For, not only do scaling studies provide accurate information on second order 
behavior, but they can be conducted in the realistic setting of lattice specific complications. What 
we have shown is that both finite lattice spacing and volume effects imply well specified correction 
terms to the scaling relations. These corrections involve both computable universal terms and 
nonuniversal terms that must be fit to a small set of measurements. Furthermore, even if one 
could conclude from lattice gauge studies that the A''^ = 2 chiral transition is second order, as 
indications show, it is still necessary to know what the infiuence of lattice specific effects are. Since 
ultimately the goal of lattice studies of finite temperature QCD is to provide accurate quantitative 
information about various properties such as transition temperatures and the fermion condensate, 
it is essential to understand the lattice specific effects. It is well known by universality that the 
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existence of a second order phase transition occurs under very general conditions, independent of 
most microscopic details. Thus for lattice studies only to establish that the transition is second 
order, although a nontrivial first step, is still not enough. The next step is to assess the accuracy 
of the results and this depends on the influence of the lattice, especially the lattice spacing. 

For this reason in many ways it is fortuitous that the chiral Nf = 2 phase transition may 
be second order, since there is a well developed theoretical framework from the renormalization 
group theory to assist in studying these delicate issues. In the present situation, where there is 
no definitive answer on the nature of the Nf = 2 transition or even whether it is a transition or 
crossover behavior, scaling theory provides a large set of consistency checks to determine if it is 
second order. 

The scenario hypothesized by |||, § assumes that the fermionic degrees solely determine the order 
parameter, neglecting any explicit influence from gluonic degrees. This is a strong assertion that 
should be placed to stringent tests. Although indications from both lattice studies and theoretical 
models further support this outcome, they leave sufficient uncertainties. On the theoretical side, 
most calculations already assume in some form that the order parameter is 0(4) symmetric and then 
proceed to compute a transition temperature by some variant of a gap equation. What is missing, 
and no doubt nontrivial, is a derivation from QCD by way of coarse graining, for example, in which 
the effective 0(4)-LGW action is derived. Turning to lattice studies, here technical uncertainties 
and results that show strong dependence on the lattice spacing jl^, ^, 21 1 obscure any definitive 



answers. 



Since QCD is dynamically much richer than most statistical mechanical systems studied on a 
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lattice, further testing would be useful in order to probe its refined properties. With little guidance 
from analytic means, the task is placed on lattice studies to first determine whether there is a 
transition or crossover effects and then obtained the detailed features of whichever. Although 
present day wisdom makes the scenario of [||, g] most plausible, as we stated above the gluonic 
degrees play no explicit role. However it may be possible that Wilson-loop-like excitations, such 
as glueballs perhaps, are sufficient in number to play a role near the reported transition. In such a 
case it is possible that the appropriate effective action which describes the loop dynamics is like the 
one described by Polyakov |22], which is also a chiral invariant theory. This could then change the 



order parameters of the theory and the transition to first order, fluctuation induced first order or 
even eliminate the transition altogether. Another possibility discussed in ||l| is that if the instanton 
density just below Tc is low enough so that r]' becomes sufficiently light, this would change the 
order parameter from 0(4) to 0(2)xO(4) and produce a first order transition. 

The above examples illustrate possible alternative interpretations, which present data can not 
rule out. To resolve between them and other such, as a first step the predicted scaling from 
renormalization group calculations should be tested against lattice simulation data. Furthermore, 
to make an accurate analysis, it would be simultaneously worthwhile to account for the effects 
of finite lattice spacing and volume. The former is especially important since it alters the single 
exponent parametrization of scaling into two. 

Apart from the fundamental issues regarding the nature of the transition, accepting that it is 
second order, there are practical instances where the lattice spacing has qualitative effects that alter 
the continuum theory much as is already known to be true at zero temperature in usual lattice 
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gauge studies. In zero temperature studies the most common example is the violation of chiral 
symmetry even for a massless theory. In similar respect, aspects of the continuum behavior of 
systems at finite temperature can be modified on a lattice. To site an example of such a case, the 
effect of finite lattice spacing introduces dimension four operator corrections to the effective action. 
This in turn prevents a tricritical point from occurring by a simple fine tuning as conjectured by 
Wilczek P]. His claim is that the effect of a strange quark is to renormalize operators of dimension 
four and higher. Thus, he asserts that by appropriately fine tuning its mass, it is possible to make 
the quartic coupling vanish so leaving the the next order, cff' interaction, as the dominant one. He 
then notes that this theory gives a tricritical point. Although lattice-model studies have been done 
on tricritical transitions in the past with no problems, the scenario pictured by Q is slightly more 
subtle. He envisions that the strange quark field induces operator insertions of degree four and 
higher so as to renormalize the original quartic coupling. This effect from the strange quark field 
of renormalizing the couplings is the same as will also occur from finite lattice spacing. This means 
in the effective theory, both effects mesh together in redefining the operator coefficients. 

Although in the continuum theory where the lattice effects are removed, there is no obvious 
problem with the so described scenario of the strange quark, if one accounts for the lattice spacing, 
quartic operator corrections from all sources will hamper the dimension six terms from controlling 
the critical behavior in a manner as simple as stated in 0. Nevertheless, if the effects of the 
unwanted quartic terms can be made small, by a small enough lattice spacing, it may be possible 
to detect the crossover towards the tricritical point. 

Another place where finite lattice spacing effects are relevant is in estimating the energy density 
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in a volume of order the correlation length near Tc- With a finite lattice spacing, the correlation 
length is modified by a factor (l+C^L{a)t^'^). Since t^'^ tends slowly to zero near Tc these corrections 
are potentially important and should be accounted for in future estimates. Also important are finite 
volume effects not only for the usual technical reasons, but also because they more closely simulate 
physical reality in one important practical applications, that of heavy ion collision. 

To summarize, in this paper we have discussed the use of the renormalization group theory of 
scaling corrections in the context of lattice gauge theory. A program based on the scaling relations 
given here would require one empirically determined nonuniversal function that depends on the 
lattice spacing, which could be obtained from fitting to a minimal set of measured thermodynamic 
quantities. Having done this, one would then have a complete set of scaling relations that include 
effects of finite lattice spacing. These relations would provide one more cross check for testing the 
consistency of simulation data. For completeness it should be clear that we have discussed only 
the leading correction from lattice spacing. Higher order corrections may also be important. This 
depends on what studies find based on this simplist extension. 
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Appendix A 

In this appendix we will show that the lattice dependences in scaling correction terms are of the 
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forms given in (3-4) to (3-9) and (3-12). In particular we will establish that these correction terms 
depend on the lattice through one nonuniversal function L{a). Our derivation closely follows those 
given in H, g]. We hope our derivation below simplifies theirs by being closer to the standard field 
theoretic language. We will need correlation functions at temperature t and magnetization M. For 
this we will have to be slightly more general than in Section 2. Following standard definitions, let 
r(^:^) denote the bare vertex function with L insertions of 0^ and N insertions of (j). This is related 
to the renormalized correlation functions of the massless theory by, 

^ R ill, ■ ■ ■ , Qn;pi, ■ ■ ■ ,pl; g, fj; ■^) = Z^2Z^''^T'^'^'^\qi, . . . ,qN;pi, ■ ■ ■ ,pl; go, J^)- (A-1) 

Following 1^ we take our renormalization conditions as, 

rg'°^('?,-g;9,/^,A)|,2=o = (A -2a) 

^rg'°)(g,-g;ff,/i,A)|,2=^2 = 1 (A -26) 

Tr ({g};ff,At,A)U.p. = g^J'^ 
^R (qi, Q2,p; g, iJ-, ■^)\s.p.' = 1 
^R {p,-p;g,tJ',-^)\p2=f,2 = 

where s.p. means in addition to q"^ = |/x^ that qi ■ qj = |/x^(4(5jj — 1). The point s.p.' is where 
q^ = jfj.'^,qi ■ qj = —\li^ so that p^ = {qi + 52)^ = /^^- We can construct correlation functions at 
finite t and M using the ones above, which were for the symmetric, massless theory, as follows, 

r(^'^)(g„p,;t,M,9,/i)= ^ 'r^+l^L^J(q^,p^.,g,^) (A - 3) 

i,j=o ^■'^■ 
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where t is related to ttiq of the bare theory by, 



2 2|/2 2\ 2 , ry , 

rriQ = m^ + {m^ — m^) = m^ + Z^2t 



(A -4) 



so that i — > corresponds to the critical point. 

Our task is to determine the scaling behavior of the correlation functions in terms of the 
temperature and relate the correlation functions to thermodynamic quantities. The renormalized 
correlation functions satisfy the renormalization group equation, 






N + M 



d 



Ha) 



L + 



d_ 

di 



)rif'^^(g,p,<7,t,M)=0 (A -5) 



where g = ^ ''u\s the dimensionless coupling constant and 



m 



dg 



-1 



1 , .91nZ«; 

2 = f3{g)^-^ 



<g) 



V{9) = P{9) 



dlnZ. 



which excludes the case L = 2, A^ = at nonzero magnetization, where an additional term is 
needed on the right hand side. This quantity will not be considered in our work. Above, ^fo is the 
dimensionless bare coupling constant, which by dimensional analysis can only depend on g and the 
ratio — . In the limit A -^ oo, go will therefore only depend on g. This fact allows one to express 
P{g) in the form given above. Further discussion of this careful point is given in chapter 8 of [^. 
To obtain the homogeneous solution of (A-5), one introduces a parameter A such that 



In A 



9W dg' 
W) 



(A - 6a) 
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with g{X = 1) = g and , 



t{X) = texp 
M{\) = Mexp 

'- R 



fiiX) = \n 



9 ^Hg') ^Pig') 



1 f9W7]{g') , 



iq,p;g,t,M) 



(A - 66) 
(A - 6c) 
(A - 6d) 
(A - 6e) 



i?^(A,5)5^(A,5)i^'^('-'^(''-'+^')/')~^^if•^^((7t7^1,t-'^;5(A),t = l,M,t7-^^^^ 



from which one gets, 



r(^'^) {q,p,g,t,M,fi) = 



(A-7) 



.{N,L) q p , 






t(A) M(A) 



V'Am' 'AV'(A^)i 



dZT'^ = l 



where 



^(A,5) =exp 
S{X,g) =exp 



1 1\ dg' 



3 U(5') i^J f3ig') 



1 r^CA) f^^' 

(^(5)-^)- 



Pig') 



Since A is arbitrary, we fix it such that A ^ corresponds to the critical region t ^ 0. For this we 



make the standard choice 



t(A) 
AV^ 



(A 



so that from (A-6c) when A ^ we have 



A=[i?(A,5)C = t 



r 



(A -9) 



30 



We want to examine the solution (A-7) for g near g* in the direction specified by the trajectory 
along A. Thus we Taylor expand (A-7) with respect to g as, 



TRig) =TRig*) + (giX) - g*)^\g=g. 



(A -10) 



From (A-6a) and using l3{g{X)) = ij{g{X) — g*) we find, 



1 



din A 



dg{\) iv{g{X) - g* 



(A -11) 



where u! = (3'{g*) is the Wegner exponent Q. From this we have. 



g{X)=g* + {g-g*)X^ . 



(A -12) 



Substituting into (A-10) we get the expansion, 



rif'^)(g,p;ff,t,M,/i) = 5^i?it-K«(-^{^-2+.)/2)-L 



^if'^)(5t,-^pt;'^;g^t = l,SMt~'^('^-2+'')/2,^ = 1) 



l + {g-g*)R--t 



f^^^'^\gt--,pt--,gM,gMt--('^-^+'^) 



(A-13) 
,/^=l) 



^(^'^)(gt^^lJtr"^5^l,^Mt 



-v{d-2+r]) 



11 = 1) 



-Co] 



All lattice dependence is contained in g. In particular this means the correlation functions have an 



overall lattice dependent amplitude factor R and a single function L{a) = g — g* that enters as a 



factor in the leading scaling correction term. 

What remains is to relate the correlation functions to the thermodynamic quantities of interest. 
The correlation length is defined as. 



.2 _ grK/V|,.^o 

p(2,0)| 

'- bare \p^=0 



r(2,0)| 



(A -14) 
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so that using (A- 13) and suppressing unnecessary indices we get 



.2 _ -^iP'r '^)l^-0^^^_^,^^^,^. 



g2p(2,0)(p^-^^g,) 



"'^"■°:ff""-' l,..o or^'^iptr.r)^ 



R-^t-^" 



g2r(2.0)(pt-^,g*)l 



^5 



Ip2=0 



(A-15) 



. , , ^^(2■o)(pf-^g*) | 



which agrees with (3-4). The susceptibility for T > T^ is, 



X = 






(A-16) 



i + ig- g^R^^r" ff„,; ^ 



which agrees with (3-5). Finahy the equation of state is, 



H{t,M,,j) s 



ar(..». _ ^arMfM^ ^ ^^ _ ^.^^^^^^aV'^J m,.,,') ,^.,,, 



9M 



5ff9M 



dgdM 



SH{U-i,Mr,g*) 



l + {g-g*)R^-t 



dgdM 

ar(0'0)(trMr,g*) 

dM 



where 



and the susceptibihties are, 



Mr = SM , 



H{tr,Mr,g*) = f 



-M/ 



m 



dH 

"" dM 



H 



S'mi{tr,Mr,g* 



l + {g- g*)R't 






""'"m"^'""'^*'''^"'^*^ 



—g^p-[lr,IVlr,g , 
g2p(0,0) 



^ + i9-g)Rt 

dM V '■' ''' " ^ 



(A - 186) 



(A - 18a) 
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Since the parameter A in our analysis was arbitrary, above we conveniently chose it to study 
scaling behavior with respect to the temperature t. In a similar way, one can study scaling with 
respect to the magnetization M in the equation of state. In particular when a; ^ 0, so that M ^ t, 
we have f{x) ~ 1 and the M-dependent scaling correction term dominates over the t-dependent 
one. To examine this we chose A such that (M(A))/((A;u)^'^^^ = 1 which implies 



SM 



d-2 + r) 



(A -19) 



Substituting into (A-13) we find, 






* Present address: Department of Physics, the Pennslyvania State University, 104 Davey Lab- 
oratory, University Park, PA 16802-6300 
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